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Abstract. New halo-independent lower bounds on the product of the dark matter- 
nucleon scattering cross section and the local dark matter density that are valid for 
annual modulations of dark matter direct detection signals are derived. They are ob¬ 
tained by making use of halo-independent bounds based on an expansion of the rate 
on the Earth’s velocity that were derived in previous works. In combination with as- 
trophysical measurements of the local energy density, an observed annual modulation 
implies a lower bound on the cross section that is independent of the velocity distribu¬ 
tion and that must be fulhlled by any particle physics model. In order to illustrate the 
power of the bounds we apply them to DAMA/LIBRA data and obtain quite strong 
results when compared to the standard halo model predictions. We also extend the 
bounds to the case of multi-target detectors. 
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1 Introduction 

In order to analyze a direct detection (DD) signal [1-3] it is typically assumed that 
the dark matter (DM) has a Maxwellian distribution of velocities truncated at the 
galactic escape velocity. This is known as the Standard Halo Model (SHM). However, 
A-body simulations indicate a more complicated DM halo [4-6]. The dependence of 
the WIMP signals on the velocity distribution has been studied in several works, see 
for instance refs. [7-13]. In particular, the assumed velocity-distribution is crucial 
for the interpretation of annual modulation signals [2, 3], with different unvirialized 
components such as tidal streams [14-17] or debris flows [6] modifying the results 
drastically [12]. Due to these large uncertainties, halo-independent methods have been 
developed to compare different DD signals [18-37]. In ref. [38] they have been extended 
in order to relate a positive signal in DD with the neutrino rate from DM annihilations 
in the Sun, while in ref. [39] a halo-independent upper limit on the scattering cross 
section was obtained by combining upper limits from DD and neutrino rates from the 
Sun. 

Recently, in ref. [40] a new framework has been established that allows to compare 
the constant rate of a DD signal with local energy density measurements, LHC limits, 
the relic abundance and indirect detection constraints (see also ref. [39] for another 
method to derive a lower bound on the cross section). In this work we show how the 
observation of an annual modulation signal in a DD detection experiment can also be 
used to place a lower bound on the product of the local DM density and the scattering 
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cross section that is independent of the DM velocity distribution, in a similar fashion 
as was done in ref. [40] for constant rates. For this purpose, we use halo-independent 
bounds on the annual modulation based on an expansion of the rate on the Earth’s 
velocity [24, 25]. In order to illustrate their use we perform a numerical analysis to 
DAM A/LIBRA data, which we will denote by DAMA in the following. As we will see, 
in combination with measurements of the local energy density we obtain quite strong 
lower bounds on the cross section, that must be fulhlled by any particle physics model 
in order for the signal to be consistent with DM. 

The paper is structured as follows. In section 2 we set the notation for DD and 
the annual modulation. In section 3 we review the case of DAMA modulation in the 
context of the SHM and we perform a £t to the signal. This section can be skipped 
by the readers familiar with the DAMA signal. Section 4 encodes the main results of 
this paper, where we derive lower bounds on the product of the local DM density and 
the scattering cross section that are independent of the DM velocity distribution and 
are valid for annually modulated signals. We do this in two ways: in section 4.1 by 
making use of the fact that constant rates need to be larger than modulations, and in 
section 4.2 by using an expansion of the rate on the Earth’s velocity [24, 25]. As an 
illustration we apply our halo-independent bounds to the DAMA signal in section 4.3. 
Finally, we give our conclusions in section 5.^ 


2 Dark matter direct detection and the annual modulation 
signal 


In this section we review the relevant expressions for DD. We focus on elastic scattering 
of DM particles y off a nucleus with mass number A depositing a nuclear recoil energy 
Bji. The differential rate for a detector with n different target nuclei is given by: 


TZ(E^,t) = E^Y.— f 
^ 


— I d'^vvfdet{v,t)^^^{v) , 


dE 


= 


R 


IrriAER 


( 2 . 1 ) 


where is the local DM mass density and is the minimal velocity of the DM particle 
required for a recoil energy Er, with p^A being the reduced mass of the DM-nucleus 
system. Ja is the mass fraction of the detector nucleus labeled by A. For single¬ 
target detectors, there is just one contribution and thus the sum over A is absent. 
fdet{,v,t) describes the DM velocity distribution in the detector rest frame, with the 
normalization J d^v fdet{v,t) = 1- The velocity distributions in the rest frames of 
the detector, the Sun and the galaxy are related by fdet{.v,t) = fsuniv + Ve{t)) = 
/ga^T + Fs -|- Veit )), where Veit) is the velocity vector of the Earth relative to the Sun 
and Vs is the velocity of the Sun relative to the galactic halo. 

We will concentrate on elastic spin-independent (SI) and spin-dependent (SD) 
contact interactions, where the differential scattering cross section daA{v)/dER scales 

^In appendix A we give a detailed derivation of the upper bound on the modulation. In appendix B 
we show the binned version of the bounds, while in appendix C we derive them for multi-target 
detectors. 
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as 1/u^, and we will take equal couplings of the DM to neutrons and protons. Then 
the SI cross section becomes^ 


daA 

dEji 


(n) 


771 ^( 731 ^^ 


FI{Er) , 


( 2 . 2 ) 


where cxsi is the total DM-proton scattering cross section at zero momentum transfer, 
li^p is the DM-proton reduced mass, and Fa{Er) is a nuclear form factor. For SD 
interactions there is no enhancement, the form factor is different, and (Tsd will 
denote the zero-momentum DM-proton scattering cross section. 

Using eqs. 2.1 and 2.2 the event rate can be written in a very compact way as 


7^(U«, t) = Y, Fa{Er, t) = Cj2U A^EHEr) t) , (2.3) 

A=1 A=1 


where we dehned the detector-independent quantity 


(2.4) 


Due to the time-dependence of the Earth’s velocity and its orientation with respect to 
the Sun’s velocity DD signals exhibit annual modulations [2, 3]. The time dependence 
in the event rate is introduced through 


+ (2.5) 

where we have dehned the constant contribution as 



and the modulated one comes from an expansion to hrst order in Ve = 29.8 km/s 
<C Vsun — 230 km/s: 

= Ap{v^)cos27i[t-E{ER)]. (2.7) 

Here An{v^) = 0.5 [6r]{v^,t*{ER)) - 6r]{v^,t*{ER) + 0.5)], where the phase t^{ER) is 
the time of year at which the event rate is maximum, and thus Ap{v^ > 0. Notice 
that observing a sign hip in the modulation amplitude would be a strong evidence that 
the signal is due to DM. However, the phase and the energy at which the hip may 
occur depends on the velocity distribution. Moreover, there could also be several sign 
hips at diherent energies. For the SHM, = 0.42 y, corresponding to June 2nd, 
and there is a sign hip at ~ 210 km/s, see for instance ref. [12]. 

Similarly, the amplitude of the modulated rate is: 

A4{Er) = - [JZ{Er, t^) — 1 Z{Er, E + 0.5)] , (2.8) 

^In our notation “A” specifies the nucleus and represents its mass number. 


3 








which is always positive. 

For a specific detector the number of DM induced events in an energy range 
[F^i, E 2 ] is given by 

"\EuE,](t) = MTCl[r,(v*,t)]]E,, (2.9) 

where M and T are the detector mass and exposure time, respectively, and we define 

^ nOG 

= dERFHER)G(i,^^„^,{En)X, ( 2 . 10 ) 

24=1 >^0 

which is an integration over energy of the weighted sum of the target nuclei for the 
quantity X{v^). Here is the detector response function describing the 

probability that a DM event with true recoil energy E^ is reconstructed in the energy 
interval including efficiencies, energy resolution, and quenching factors. 

3 The DAMA annual modulation signal within the SHM 

In section 4 we will derive halo-independent bounds valid for annual modulation signals, 
and in order to illustrate their use, we will apply them to the modulation observed 
by DAMA [41]. In the following we review and analyze the interpretation of this 
modulation within the SHM, and therefore this section can be skipped by those readers 
familiar with this signal. 

The DAMA experiment reports an annual modulation of the rate in their Nal 
scintillator detector, with a period of one year and a maximum around June 2nd. With 
an exposure of 1.33tony the modulation is claimed to be compatible with DM at a 
very high statistical significance (9.3 (t) [41]. The amplitude of the annual modulation 
observed in the [2, 6] keVee range is A4[2,6] = (0.0112 ± 0.0012) counts/keVee/kg/day 
[41], while above this range the data is consistent with no annual modulation. As 
it is well known, the signal can be explained under standard assumptions regarding 
the DM halo (SHM) and WIMP interactions. However, the DM interpretation of the 
modulation is strongly disfavored by other experiments independently of the halo, both 
for elastic SI and SD interactions [25], and for inelastic interactions [28]. Despite this, 
we will use the DAMA signal in the following to illustrate the use of the lower bounds 
derived in section 4. Notice that if DAMA were not due to DM, current limits (for 
instance by LUX [42]) imply that observing annual modulations will be very hard in 
the future (see for instance ref. [43]). 

First let us discuss which is the dark mass range consistent with the data for the 
SHM. There is no observation of a sign flip and all the observed modulations in the 
[2, 6] keVee range are positive, which means that for the SHM they are above the sign 
flip, implying that Vm ^ 210km/s. Therefore, in order for the modulation of the first 
bin to be positive, if the scattering is on sodium, from eq. (2.1) the DM mass needs to 
be smaller than ~ 30 GeV, while if the scattering is on iodine, it needs to be smaller 
than ~ 90 GeV. In addition, there are also lower bounds on the DM masses in order 
for the DM particles to be energetic enough to deposit recoil energies above threshold. 
This means that if the scattering is on iodine, DM masses below ~ 20 GeV cannot 
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Figure 1. Minimum velocity VmiEji), eq. (2.1), versus the DM mass shown for sodium (red) 
and iodine (blue) for recoil energies relevant for the DAM A experiment: Eji = 2, 4,6keVee 
(from bottom to top as solid, dotted and dashed curves). It is shown as dotted black the Vm 
below which the modulation is negative for the SHM, and as dashed black the typical escape 
velocity in the detector rest frame. The shaded area encoded between the dotted vertical red 
(blue) lines is the allowed mass range within the SHM for scattering on sodium (iodine). 


produce observable recoils. On the contrary, for > 20 GeV and SI interactions, 
iodine typically dominates the rates due to the enhancement.^ Furthermore, by 
demanding that there are observable recoils in the 6 keVee range, we obtain that 
> 8 (30) GeV for scattering on Na (I). 

This is illustrated in figure 1, where we plot the minimum velocity VmiEn), 
eq. (2.1), versus the DM mass for sodium (red) and iodine (blue) for recoil ener¬ 
gies relevant for the DAMA experiment: = 2, 4,6keVee (from bottom to top as 

solid, dotted and dashed curves). We also show as dotted black the minimum velocity 
below which the modulation flips sign for the SHM, and as a dashed black line the 
typical escape velocity in the detector rest frame. The shaded area encoded between 
the dotted vertical red (blue) lines shows the allowed mass range 8 < (GeV) < 30 
(30 < (GeV) < 90) for scattering on sodium (iodine) within the SHM. For more 
complicated haloes with streams or debris flows the phase of the modulation will be 
different (see for instance ref. [12]), and so will the allowed mass range. 

In order to analyze more quantitatively DAMA’s modulation within the SHM, 
we perform a Sf fo the data in the [2,6] keVee energy range, both for SD and SI 
elastic interactions. For our numerical analysis we take equal couplings to protons 
and neutrons, both for SI and SD,"^ and quenching factors gNa = 0.3 for sodium and 

^For this reason lowering the threshold of the DAMA experiment would in principle allow to 
distinguish between the two SI solutions, the ~ 10 GeV DM mass, for which a lower threshold would 
imply an enhanced rate due to scatterings on iodine, and the ~ 80 GeV one, which would show a 
phase flip in the lowest energy bin, see ref. [44]. 

^Regarding SD, both ^^Na and ^^^I are dominated by the spin of the protons. Therefore, suppressed 
couplings to neutrons help to reduce the inconsistency with other null-result experiments that have 
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I 

Na 


(GeV) 

CTSI/SD (cni^) 

xLn/dof 

(GeV) 

(Tsi/sD (cm^) 

xLn/dof 

SI 

79.4 

1.1 ■ 10-^^ 

7.7/6 

12.6 

1.8 ■ lO-^^* 

8.3/6 

SD 

63.1 

5.0 ■ 10-37 

7.9/6 

12.6 

6.3 ■ 10-37 

8.7/6 


Table 1. Results of the DAMA fit to the SHM, done for SI and SD interactions, assuming 
equal couplings to protons and neutrons. We use Uesc = 550 km/s, = 0.4 GeV cm“^ and 
the quenching factors ^Na = 0.3 for sodium and qi = 0.09 for iodine. 


qi = 0.09 for iodine. For the SI form factor we use the Helm parameterization, F{Er) = 
3e-9 « /2[sin(gr) — gr cos(q'r)]/(gr)^, with = 2mAEji and s = 1 fm, r = \/E? — 5s^, 
R = gp we take the structure functions from ref. [46]. For the 

SHM we use a Maxwellian velocity distribution with mean velocity v = 220 km/s 
truncated at the galactic escape velocity Uesc = 550 km/s, and a local energy density 
= 0.4 GeV cm“^, see ref. [47] for a recent review. 

We use the latest results from hgure 8 of ref. [41]. As expected from the literature, 
there are two minima: for the large DM mass solution the scattering is predominantly 
on iodine, while for the low DM mass one it is on sodium. The best-£t values are shown 
in table 1. In all cases we obtain a good £t, with Xmin/dof < 1.5, where dof stands for 
the number of degrees of freedom, which is 6 (8 data points minus 2 parameters). Notice 
that for SI interactions, for the large DM mass solution the scattering is predominantly 
on iodine and therefore the cross section is much lower than in the case of low DM 
masses (scattering on sodium), due to the larger enhancement of the former. For 
SD both solutions imply roughly the same cross sections.^ Our results agree with those 
present in the literature, see for instance refs. [45, 48, 49]. 

In hgure 2, we show the allowed 90% CL parameter space for the SHM as solid 
blue lines, together with the = 2.3, 5.99, 9.21 contours (in blue, green, light blue 
respectively) corresponding in two dimensions to a CL of 68.27, 95, 99% respectively. 
This is done for SI (top) and SD (bottom), and we show the low (high) mass solution 
in the left (right), which correspond to scattering on Na (predominantly on I). We also 
plot the best-£t points with blue marks. 


4 Halo-independent lower bounds on p^^dsi/SD for annual mod¬ 
ulations 

In this section we derive halo-independent bounds valid for annual modulation signals. 
It encodes the most important results of this work. 


nuclei dominated by the spin of the neutrons, such as ^^Ge, ^^®Xe and ^^^Xe, see for instance ref. [45]. 

^This behaviour can be clearly seen in figure 5 of appendix C, where we show in the same plot 
both solutions in blue (dotted for scattering on Na, dashed on I, and solid on both), for SI (top) and 
SD (bottom), together with the confidence level regions. 
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Figure 2. SHM best-fit regions for the DAMA data for SI (top) and SD (bottom), and 
scattering on Na (left) and I (right). The = 2.3,5.99,9.21 contours corresponding in 
two dimensions to a CL of 68.27, 95, 99% are shown in blue, green and light blue, respectively, 
together with the best-fit points, depicted with blue marks. 

4.1 Using that the constant rate is larger than the annual modulation 

An upper bound on the halo integral r]{vm)^ can be easily obtained using the normal¬ 
ization condition of the velocity distribution as well as the definition of eq. (2.6), 

as was done in refs. [32, 36, 40]: 


poo poo 

1 = / d^vfdet{v) = / rj{v)dv 

Jo Jo 

(4.1) 

>fj{vi)vi+ / dvfj{v) 

J Vi 

(4.2) 

> fj{vi)vi , 

(4.3) 


where from the hrst to the second line we used that below the threshold vi the minimum 
is obtained for constant fj{y) = 17 (^ 1 ), as it is a monotonously decreasing function, and 
in the last line we dropped the second term. 

From these expressions an upper bound on the constant rate can be derived [40]. 
In order to have a bound valid for annual modulation signals, one can simply use that 

®To alleviate notation, in the following we will drop the A superscipt for the minimum velocity, 
v:^, see eq. (2.1). We will show it explicitly when it is crucial to illustrate its dependence on the 
nucleus. 
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the constant rate is larger than the amplitude of the modulation: 


^ -^ 77 (^ 771 ) • 


(4.4) 


If the spectrum is not measured, one can derive a bound in terms of the measured 
number of modulated events in an energy interval [Ei, E- 2 \ by combining eq. (2.9) with 
eqs. (4.3) and (4.4). Using the dehnition of C, eq. (2.4), one gets: 


Px'^si > 


MT iii/t-^yil; 


AT-mod 

^^[Ei,E2] 


“Events bound”, 


(4.5) 


where |ll/<]|®( is dehned in eq. (2.10). Notice that this bound can be used for multi¬ 
target experiments. 

If the spectrum is measured, one can obtain a stronger bound using eq. (4.2) 
instead of eq. (4.3). In this case, for a single-target experiment with perfect energy 
resolution, for hxed DM mass, a measurement of the spectrum A4 {Er) allows a deter¬ 
mination of the halo integral via eqs. (2.3) and (2.8): 


-^ri (ym) 


M{Er) 
CA^EHEr) • 


Then, in combination with eq. (4.2), the following bound can be obtained: 


(4.6) 


Px^si > 


2"ixhxP 


M{E,) 

"’"niEi) 




fvi 


Fl 


{Er)J 


“Spectrum bound”. (4.7) 


If a DD experiment observes an annual modulation spectrum, eq. (4.7) provides a lower 
bound on the product Px'^si which is independent of the local DM velocity distribution. 
In ref. [40] bounds valid for constant rates analogous to those of eqs. (4.5) and (4.7) 
were derived. Notice that if a sign flip is observed, a better approach is to apply 
the bound for the region below and above the sign flip separately, in order to obtain 
stronger bounds. 

If the detector has different elements, in order to use eq. (4.7) one needs to assume 
that one of them gives the dominant contribution. The bound will be stronger in 
general for the lowest mass nuclei in the case of SI, due to the A^ suppression, but so 
will the SHM prediction, and therefore we expect the ratio of the bound and the SHM 
to remain approximately the same for different nuclei. One also needs to take into 
account that, although for SI interactions the heaviest element typically dominates the 
rate due to the A^ enhancement, for low DM masses the scattering on the heavy nuclei 
may not be possible, as there are no DM particles with speeds larger than the escape 
velocity in the detector rest frame, ~ 750kms“^. We will illustrate these features for 
the DAMA data in section 4.3. In appendix C we generalize eq. (4.7) to the case of 
multi-target detectors. 







4.2 Using upper bounds on the annual modulation 

One a can derive stronger tests for the case that the spectrum is measured by using 
the results of refs. [24, 25], where upper bounds on the annual modulation signal in 
terms of the constant rate were obtained by doing a hrst order expansion of the rate 
on the Earth’s velocity. 

Lets consider a spectral measurement of both the modulation A4.{Er) and the 
constant rate TZ{Er) in the energy range [^ 1 ,^ 2 ], which for a given DM mass can 
be related to a velocity interval [^ 1 ,^ 2 ] via eq. (2.1). By doing an expansion on the 
Earth’s velocity and some mild assumptions about the halo (see below) one can derive 
the following upper bounds on the modulation amplitude [24, 25]: 



“General halo”, (4.8) 


“Symmetric halo”. (4.9) 


For the “General halo” it is assumed that the velocity distribution is constant in time 
on the scale of 1 year (i.e., that all the time dependence of the rate comes from the 
velocity of the Earth) and that it is constant in space on the scale of the size of the 
Sun-Earth distance. Notice that halo substructures like streams with velocities larger 
than Vf. are covered by the bounds. The “Symmetric halo” further assumes that there 
is a preferred direction for the DM flow, specified by a, the angle between the DM 
flow and the orthogonal to the ecliptic. The most conservative bound is obtained for 
sin a = 1, which corresponds to a DM stream parallel to the ecliptic. However, it may 
happen that the DM velocity is aligned with the motion of the Sun, and therefore 
sin a ~ 0.5. This is the case of isotropic velocity distributions and, up to a small 
correction due to the peculiar velocity of the Sun, it also holds for triaxial halos or 
a dark-disc. For the “General halo” the bound is independent of the phase, which is 
free. However, for the “Symmetric halo” the phase is independent of Vm (and therefore 
independent of Er), and fixed. The SHM is included in this last case, with sin a ~ 0.5 
and the phase being equal to June 2nd. As for the case of the “Spectrum bound”, 
eq. (4.7), if a sign flip in the modulation is observed in the velocity range considered, 
one may apply the bound for the region below and above the sign flip separately, in 
order to get stronger bounds. More details regarding the different assumptions of the 
bounds can be found in the original ref. [24] (and also in refs. [25, 28]). 

From the above discussion it is clear that in order to decide whether to apply the 
“General halo” or the “Symmetric halo” one can use the information encoded in the 
phase. First of all, one needs to do a fit keeping as free parameters the modulation 
amplitude, the period and the phase. If data disfavours a constant phase across the 
different energy bins, this invalidates the assumption of a “Symmetric halo” and one 
is forced to use the “General halo”. Regarding the polar angle a, unfortunately, the 
phase of the modulation does not carry information about it, and different values of a 
will give the same phase. It only depends on the azimuthal angle within the ecliptic, 0, 
on which however the bounds do not depend. This is easy to visualize if one imagines 
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a stream parallel to the ecliptic, with a = /2. In this case it is clear that the phase 

is completely determined by the angle 0 at which the Earth’s velocity and the stream 
direction are aligned. Notice that the amplitude of the modulation is proportional to 
sin a, but, unfortunately, its value cannot be disentangled without specifying the halo. 

In order to check the consistency of a DD signal one can apply these bounds if 
both the modulation and the constant rate have been measured [24]. In some cases one 
may not be able to separate the background from the DM contribution in the measured 
constant rate, as is the case of DAMA, but it is clear that if the inequalities are violated 
for an unknown constant background plus DM rate, they will also be for just the DM 
component, and the conclusion (the exclusion of the DM as an interpretation of the 
modulation) will still hold. However, if the bounds are fulhlled for some unknown 
constant background plus DM rate, they may well be violated for just the DM rate. In 
these cases, in order to check if the modulation observed in one experiment is consistent, 
one can compare it with upper limits on the rate from a different experiment [25, 28]. 
Interestingly, as we will show now, one can obtain upper bounds that do not depend 
in any way on the constant rate nor on its upper limits. In combination with other 
measurements, these set halo-independent constraints on the DM models [40]. 

By applying eq. (4.2) in the “General halo”, eq. (4.8), we can obtain an upper 
bound on the modulation (see appendix A for the derivation): 


/ 2 1 

dvAJv) < Ve { - 

V2 


'VI 


(4.10) 


Similarly, for the “Symmetric halo”, eq. (4.9), we obtain: 


rV2 

dvAn{v) < sin a Ve — . 


>Vi 


(4.11) 


Both bounds are written completely in terms of velocities, with no factors containing 
the constant rate. Therefore, even if just the annual modulations are measured, one 
can get an upper bound. This is crucial, since, as we argued, some experiments are 
not able to disentangle constant DM signals from constant backgrounds. 

If we now consider a spectral measurement of the modulation such that we can 
extract via eq. (4.6), then, using eq. (4.10), we can get an upper bound on the 
integrated modulation over the velocity range: 



M(v) 

A^Fi(En) 


<CVe 




(4.12) 


from which we can derive a lower bound in the product of the energy density and the 
cross-section which does not depend on /(u), and is valid for very generic haloes: 


Px^si > 


2?»xhp ^ 
A2 Ve 



-V 

V2) Jvi Fa{Fr) 


“General bound”. (4.13) 
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An analogous bound can be obtained for symmetric haloes, eq. (4.11):^ 


Px'^si > 



1 \-i 1 


sm a Ve 


1)1 


dv 


M{v) 

FI(Er) 


Symmetric bound”. (4.14) 


If a DD experiment observes an annual modulation, eq. (4.5) (“Events bound”), 
eq. (4.7) (“Spectrum bound”), eq. (4.13) (“General bound”) and eq. (4.14) (“Sym¬ 
metric bound”) provide lower bounds on the product which are independent of 

the local DM velocity distribution. The first two just use that modulations are smaller 
than constant rates, while the last two are based on an expansion of the rate on the 
Earth’s velocity, and are clearly a factor ~ v/v^ (with v > Vm) stronger than the 
former. These are the main results of this paper. Notice if the detector has different 
nuclei, in order to use eqs. (4.7), (4.13) and (4.14), one needs to assume that one of 
them gives the dominant contribution. In appendix C we generalize these bounds to 
the case of multi-target detectors. 

4.3 Applying the halo-independent bonnds: the case of DAMA 

Now we will apply the halo-independent bounds derived in the previous section to the 
annual modulation observed by DAMA [41] (see section 3 for a review and analysis of 
DAMA in the context of the SHM). As discussed in section 4.2, in order to know if one 
can apply the “Symmetric bound” or only the “General bound” one should perform 
a fit to the modulation data leaving the amplitude, the period and the phase free in 
order to determine whether the phase is constant or not across the energy range of 
the signal. For DAMA these fits have been performed [41], obtaining that the phase 
is compatible with being constant across different energy bins within 3 ct, and equal to 
June 2nd (see table 4 of ref. [41]). 

It should be emphasized that the DAMA modulation is compatible with the SHM 
(and thus with sin a ~ 0.5), see section 3, but it can also be caused for instance, as 
discussed in section 4.2, by a stream with a different polar angle, sin a ^ 0.5, but with 
the same azimuthal angle (j) as in the SHM (where it is given by the Sun’s velocity) such 
that the phase is June 2nd, and with a stream speed such that the hnal modulation 
amplitude has the same value as the one observed by DAMA. Therefore, for DAMA 
one can use the “Symmetric halo” for different values of a. In this case, the most 
conservative option is sin a = 1, and if one is interested in analysing a DM flow aligned 
with the velocity of the Sun one can use sin a ~ 0.5. In order to illustrate their 
different strengths we will also show the “General halo”, which does not require a 
constant phase.® 

^We keep it in this form of instead of writing vi for easier comparison with the binned 

version, see eq. (B.5) in appendix B. 

®Notice that the differences in the modulation amplitudes and their errors in the case of fixed (free) 
phase, c.f. table 3 and 4 of ref. [41], are negligible, and therefore we can use the same data when 
applying both bounds. However, in general this may not be the case, see for instance the CoGENT 
analysis of ref. [24], which uses the extracted modulation amplitudes for the different assumptions of 


ref. [50]. 
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In order to apply the bounds to real data we need binned expressions, which 
are given explicitly in appendix B. To extract A^{vm) from the observed modulation 
amplitude M.{Er) via eq. (4.6), and therefore to apply our bounds, one needs to assume 
that the scattering occurs on a particular nucleus. We will assume in the following 
that the scattering is dominated by either sodium or iodine, which, as we explained 
before, is a very good approximation, given the large mass hierarchy between them. 
In appendix C we generalize the bounds to multi-target experiments and apply them 
to both nuclei at the same time. 

In hgure 3 we show the 90% CL lower bounds on the scattering cross section, 
together with the SHM preferred regions, already shown in hgure 2. For dehniteness, 
we use the mass range allowed for the SHM (plotted in blue), see section 3 for details. 
From bottom to top, we plot the “Spectrum bound”, eq. (4.7), which just uses that 
Ari < r], as solid red; we also show the bounds that use the expansion on the Earth’s 
velocity: the “General bound”, eq. (4.13), as dashed red, and the “Symmetric bound”, 
eq. (4.14), as dotted (dotted-dashed) red for sina = 1 (sina = 0.5, i.e., for the DM 
how aligned with the velocity of the Sun). As expected, in all cases the latter bounds 
are much stronger than the “Spectrum bound” (solid red), due the v/ve enhancement 
(c.f. eqs. (4.13) and (4.14) to eq. (4.7)). 

One can use these lower bounds on the cross section to compare with LHC limits, 
with relic abundance constraints or with indirect detection limits [40]. This has to 
be done within the context of a particular particle physics model and can be used to 
check its validity. Another useful and general comparison can be done by combining 
the lower bound on for a given DM mass with local density measurements [40]. 
This combination provides a lower bound on the cross section of any given particle 
physics model. In the analysis, we use values for in the range between 0.2 and 
0.6 GeV/cm^ (see ref. [47] for a recent review), and we hx the DM mass to be equal 
to the best-ht values obtained in the fif; see table 1. In this case one would need to 
obtain the DM mass by other means, for instance from LHC measurements, from an 
indirect detection signal (i.e., a gamma ray line) or from more than one DD signal, see 
for instance refs. [19, 37]. Notice that results for different DM masses can be obtained 
by a simple rescaling of hgure 3. 

This is shown in hgure 4, where we plot the local energy density versus the 
scattering cross section. We show in blue the 90% CL SHM predictions and in red the 
90% CL lower bounds, for SI (top) and SD (bottom), and for scattering on Na (left) 
and I (right). We plot from left to right the “Spectrum bound” of eq. (4.7) (solid 
red), the “General bound” of eq. (4.13) (dashed red) and the “Symmetric bound” of 
eq. (4.14) (dotted red for sin a = 1, dotted-dashed red for sin a = 0.5). As expected, 
the last three give much stronger constraints than the former. 

For the very conservative “General bound” (dashed red) we obtain that SI cross 
sections cxsi < 8 ■ 10“^^ (3 ■ lO”"^^) cm^ are disfavoured for = 12 GeV, scattering on 
Na = 79 GeV, scattering on I). For SD usd ^ 3-10“^® (2 x 10“^®) cm^, for = 12 
GeV, scattering on Na (m^ = 63 GeV, scattering on I). A bit stronger constraints are 
obtained for the “Symmetric bound”, eq. (4.14) (dotted red), for sin a = 1 (if the DM 
how is aligned with the velocity of the Sun, i.e., sin a = 0.5, the limits, shown with a 
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Figure 3. Results for the DAMA data for SI (top) and SD (bottom), and scattering on Na 
(left) and I (right). The SHM allowed parameter space at the 90% CL is encoded between 
the solid blue lines. The 90% CL lower bounds are shown in red from bottom to top: the 
“Spectrum bound” of eq. (4.7) (solid red), the “General bound” of eq. (4.13) (dashed red) 
and the “Symmetric bound” of eq. (4.14) (dotted red for sin a = 1, dotted-dashed red for 
sin a = 0.5. Notice that the latter is a factor of 2 stronger than the former.). The SHM 
Ax^ = 2.3, 5.99, 9.21 contours corresponding in two dimensions to a CL of 68.27, 95, 99% are 
shown in blue, green and light blue, respectively, together with the best-fit points, depicted 
with blue marks. 

dotted-dashed red curve, are a factor of 2 stronger than in this last case). Compared 
to the SHM, these are roughly between ~ 0.5 and ~ 1.5 orders of magnitude weaker. 

Notice that these cross sections are excluded independently of the unknown ve¬ 
locity distribution. In addition, the limits are conservative in the sense that they still 
hold if the particle that gives rise to the DD signal constitutes only part of the total 
DM, since then the lower bound on the total local DM density can only be larger. 

5 Conclusions 

In this paper we have derived new halo-independent lower bounds on the product of 
the DM-nucleon scattering cross section and the local DM density that are valid for 
annually-modulated signals, extending the work done in ref. [40] for constant rates. 
First, we have used that the amplitude of the annual modulation should be smaller 
than the rate. This allows to obtain the “Events bound”, eq. (4.5), if only the number 
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Figure 4. Lower bounds on the local energy density applied to DAMA modulation SI (top) 
and SD (bottom), and for scattering on Na (left) and I (right). The DM masses in each 
case are taken to be close to the best-fit points. We show the SHM prediction (solid blue), 
together with the bounds in red, from left to right: the “Spectrum bound” of eq. (4.7) (solid 
red), the “General bound” ofeq. (4.13) (dashed red) and the “Symmetric bound” ofeq. (4.14) 
(dotted red for sin a = 1, dotted-dashed red for sin a = 0.5. Notice that the latter is a factor 
of 2 stronger than the former.). All are shown at 90% CL. 


of events is measured, and the “Spectrum bound”, eq. (4.7), if the spectrum is observed. 
Then we have derived new tests by making use of halo-independent bounds obtained 
in refs. [24, 25] by doing an expansion of the rate on the Earth’s velocity. These are the 
main results of the paper: the “General bound”, eq. (4.13), valid for general haloes, 
and the “Symmetric bound”, eq. (4.14), which assumes some preferred direction for 
the DM velocity. These bounds can also be extended to multi-target detectors, see 
appendix C. 

In order to illustrate their use, we have applied them to DAMA data. In combi¬ 
nation with local energy measurements, we are able to exclude cross sections that are 
roughly between one and one and a half orders of magnitude smaller than the SHM 
ones, see hgure 4, but with the reward of being independent of the unknown velocity 
distribution. It is important to stress that these limits must be fulhlled by any parti¬ 
cle physics model that gives rise to elastic SI/SD interactions in order for the DAMA 
signal to be consistent with DM. 
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In the spirit of being as conservative as possible, we would like to encourage 
the community to show, for positive DM direct detection signals, in addition to the 
typical SHM best-£t regions, the halo-independent lower bounds derived in this work, 
for instance the “Symmetric bound” of eq. (4.14) (dotted red) for the cases where the 
phase is constant, or the “General bound” of eq. (4.13) (dashed red) when the phase 
varies, in a similar fashion as was done in hgures 3 and 5 for the DAMA signal. 

Acknowledgements: I am very grateful to Thomas Schwetz for fruitful discus¬ 
sions regarding not only this work but also our previous related projects. I would also 
like to thank Mattias Blennow and Stefan Vogl for useful discussions on our recent 
related project. 


A Derivation of the upper bounds on the modulation 


We derive here eq. (4.10) valid for the “General halo”. By applying eq. (4.2) in the 
“General halo”, eq. (4.8), we get: 






r‘V2 


CV2 


1 


f«2 


/ dvAJv) <Ve[ -/ dvr](y)+ / dv— — / dv— / dv'rjiv') . (A.l) 

Jv^ \vi vi ) 

Integrating the third term, and by parts the last term, we obtain for eq. (A.l): 


^412 / 2 

dvArj{v) < Ve 


2 f'^2 r'^2 

dvriiv) + / dv-fj(v) 

V 


'VI \^1 ^2 J vi J VI 

Now we can use that the last term of eq. (A.2) obeys 

r'02 /■«2 

/ dv-fj(v) < — / dvriiv), 

Jvi V Vi 

and inserting eq. (A.3) back into eq. (A.2), we get: 

rv2 / 2 1 1 

/ dvAJv) <Ve{ -/ dvfjiv) 

Jvi \Vl V2 Vi 


(A.2) 


(A.3) 


(A.4) 


The last term of eq. (A.4) depends on r]{v) and we can conservatively drop it. We 
hnally obtain: 

f\A,{v)<vJ---) , (A.5) 

Jvi \Vl V 2 J 

which is an upper bound on the modulation amplitude expressed solely in terms of 
velocities. In the end the result is the same as if we had used in eq. (4.8) the weaker 
bound of eq. (4.3) instead of the stronger one of eq. (4.2). An analogous derivation 
can be performed for eq. (4.11). 

Notice that strictly speaking all these inequalities should read “less than” (<), 
instead of “less or equal than” (<), because: 1) eq. (4.2) (and similarly eq. (4.3)) is 
very conservative, assuming just that rj is a non-increasing function of velocity and. 
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2) eq. (4.8), and therefore eq. (4.10) (and similarly eq. (4.9) and (4.11) for symmetric 
haloes), are based on a hrst order expansion on the Earth’s velocity. Moreover, in the 
last step we dropped a positive term.® However, we decide to keep “less or eqnal than” 
in all ineqnalities to make explicit that for very special (and probably nnrealistic) 
haloes the ineqnalities conld be (at least approximately) satnrated, see ref. [24] for 
examples. In practice this means that any modnlation signal that even satnrates the 
bonnds is very nnlikely have a DM origin. 


B Binned versions of the bounds 


In order to apply the bounds to real data, we need binned expressions, which necessarily 
involve some approximations. We will discuss the validity of the binning prescriptions 
here used at the end of the section. Let us dehne the bin average of a quantity X{E) 
as 

1 rEa 

W. = ^/ dEX(E), (B.l) 

1 J Eil 


where En and Ei 2 are the boundaries of bin i and AEi = Ei 2 — En. 

We denote by At* = {M.)i the observed modulation in bin i, where f = 1,..., N. 
For a hxed DM mass, can be obtained from the amplitude of the modulation in bin 


i, i.e.. Mi, by: 




Miq 

CA^FDiSA' 


(B.2) 


where g is a possible quenching factor which we assume to be energy-independent (as is 
the case of DAMA), n* = Vm{Ei) is the minimum velocity centered in the corresponding 
energy bin using eq. (2.1), {E\)i is the nucleus form factor averaged over the bin width 
and Ja is the nucleus mass fraction if different elements or isotopes contribute to the 
rate. Similarly, = fj{vi) can be expressed in terms of the constant rate TZi. 

Now we have to perform a bin average of the bounds. This involves quantities 
like {giy) M{E^)i, which we replace by {g{y))i Mi, which is a good approximation 
whenever g{y) does not vary drastically within a bin. The integral is approximated by 
a sum over bins. Using eq. (B.2), the “Spectrum bound”, eq. (4.7), becomes: 


Px^si > 


A2 


fA{Fl)_ 


N 


M, + Y. 


qAvi 

h f-i (P'S 


M, 


where we denote by Avi the width in velocity space of bin i. 
Similarly, for the “General bound”, eq. (4.13), we get: 


(B.3) 


Px^si > 


2m^g,p 1 

M Ve 



qAvi 
Ia {Fl). 


Mq. 


(B.4) 


^Unless vi = V 2 , which is of course not an interesting case as then there is no range in velocity in 
which a modulation has been observed. 
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(B.5) 


And finally the “Symmetric bound”, eq. (4.14), reads in its binned version: 


Px'^si > 


A2 



-1 1 
1 sin a Ve 


E 


qAvi 
Ia {Flh 


M,. 


The inevitable inaccuracies involved in the binning procedure are small whenever the 
quantities within a bin do not show drastic changes. By using different prescriptions 
for the bounds, for instance l/ui, vs ( 1 /^) 1 , vs l/(n)i, we have checked that inaccu¬ 
racies due to the binning procedure are negligible. Notice also that the averaging of 
the velocity (v) can be performed analytically, as the relation with recoil energy Eji 
is known, c.f., eq. (2.1). In addition we need to extract information from the mea¬ 
sured modulation amplitude in a particular bin (of size 0.5 keVee for DAMA), which 
necessarily implies unfolding form factors and the detector response, which involves 
energy resolution, efficiencies and quenching factors. We have also checked that by 
including energy resolution in a similar fashion as in ref. [44] (see also refs. [51, 52]) the 
quantities involved in our bounds do not vary much, and therefore energy resolution 
is not expected to change the results signihcantly. 


C Bounds for multi-target detectors 


With the exception of the “Events bound”, eq. (4.5), in order to apply the rest of the 
bounds here derived to multi-target experiments one needs to assume that scattering is 
dominated by a particular nucleus. Let us now generalize the bounds to experiments 
with different n elements, labeled by A. For a similar procedure for multi-target 
bounds see ref. [24] . In this case the modulation amplitude in an energy bin i receives 
contributions from each element: Ait = We use the notation vf = 

see eq. (2.1). Summing eq. (4.2) in its binned version for all the n elements we get: 


n n N 

A=1 A=1 i=l 


^ Q I g Avf —A 

^ CA^ED.U ^ CA^EDJ^ 



mm 

A 


UhfDJa ' ^ 


N 

E 

i=l 


mm 

A 


q Avf \ 
AHFDJ^J 



(C.l) 

(C.2) 

(C.3) 


where TZi = total rate in bin i, and minA(Xj) means that we need 

to take the minimum X between all the nuclei present in the target, for each bin i. 
From the hrst to the second line we used the equivalent of eq. (B.2) for (in terms of 
constant rates), and from the second to the third line we changed the order of the sums 
and used the fact that C is detector independent, see eq. (2.4). Using that constant 
rates are larger than modulations, TZi > Mi, and the dehnition of C, eq. (2.4), we get: 


Px<^si > 


2m^qi 
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XP 


mm 


/ g(n^)] 
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x<i+ E 


mm 


/ q Avl 
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which is the “Spectrum bound” valid for multi-target experiments. 

Now we can do a similar thing for the bounds based on an expansion on the Earth’s 
velocity, the “General bound”, eq. (4.13), and the “Symmetric bound”, eq. (4.14). 
Summing eq. (4.10) in its binned version for all the n elements we get: 


1 

C 


n N 


EE 


qAvf 

aHfdja 


Mf<Ve 




By changing the order of the sums, we finally obtain: 


(C.5) 


Px^si > 


Ve 





q Avf \ 

A^inhfAj 


Mi , 


(C.6) 

which is the “General bound”, eq. (4.13), for multi-target experiments. Doing the 
same procedure in eq. (4.11), we get for the “Symmetric bound”, eq. (4.14): 


Px^si > 


sin a Ve 




q Avf \ 


Mi , 


(G.7) 


which is now valid for multi-target experiments. 

It is easy to check how these multi-target bounds reduce to the single-target ones 
for n = 1, as it has to be. Notice however that the bounds here derived are weak due 
to the crude approximations needed in going from eq. (G.2) to eq. (G.3). For instance, 
for two isotopes of the same element, the multi-target bounds should be as strong as 
the mono-target ones, but due to the simplifications made they are clearly weaker by 
more than a factor of 2. 

In figure 5 we show in black from bottom to top the “Spectrum bound” of eq. (G.4) 
(solid), the “General bound” of eq. (G.6) (dashed) and the “Symmetric bound” of 
eq. (G.7) (dotted for sin a = 1, dotted-dashed for sin a = 0.5) for DAMA data for 
SI (top) and SD (bottom), using both sodium and iodine. The “Events bound”, 
eq. (4.5), which is valid for multi-target, is plotted as solid purple. We also show the 
SHM prediction in blue (dotted for scattering on Na, dashed on I, and solid on both). 
The SHM Ax^ = 2.3, 5.99, 9.21 contours (already shown in figure 2) corresponding 
in two dimensions to a GL of 68.27, 95, 99% are shown in blue, green and light blue, 
respectively, together with the best-fit values, depicted with blue marks. As expected, 
in all cases the multi-target bounds are weaker than the ones obtained for the single¬ 
target bounds, cf. figure 3. In particular, for SI interactions in the low DM mass 
region (m^ < 20 GeV) the difference between the multi-target bounds and the SHM 
is greater than two orders of magnitude, as the former ones are suppressed by the 
iodine factor with respect to the latter, which are dominated by scatterings on 
sodium. Notice that the “Events bound”, eq. (4.5) (solid purple), is stronger than the 
“Spectrum bound” for multi-target detectors (solid black) due to the 1/n suppression 
of the latter (n = 2 in this case), see eq. (G.4). 
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Figure 5. Bounds (in black) applied to DAMA modulation for SI (top) and SD (bottom), 
shown at 90% CL. The 90% CL SHM range is plotted in blue (dotted for scattering on Na, 
dashed on I, and solid on both) together with the “Spectrum bound” of eq. (C.4) (solid), 
the “General bound” of eq. (C.6) (dashed) and the “Symmetric bound” of eq. (C.7) (dotted 
for sin a = 1, dotted-dashed for sin a = 0.5), using both Na and 1. The “Events bound” 
of eq. (4.5), which is valid for multi-target, is plotted as solid purple. The SHM = 
2.3, 5.99, 9.21 contours corresponding in two dimensions to a CL of 68.27, 95, 99% are shown 
in blue, green and light blue, respectively, together with the best-fit points, depicted with 
blue marks. 
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